Study of QGP signatures with the 
ALICE events 



K + K~ signal in Pb-Pb 



O 

o 

c 

a 

in 



A. Akindinov 1 , A. Alici 2 ' 3 , P. Antonioli 3 , S. Arcelli 3 ' 4 , M. Basile 2 ' 3 , G. Cara Romeo 3 , M. Chumakov 1 , L. Cifarelli 2 ' 3 ^ 4 , 
F. Cindolo 3 , A. De Caro 5 , D. De Gruttola 5 , S. De Pasquale 5 , A. Di Bartolomeo 5 , M. Fusco Girard 5 , Yu. Grishuk 1 , 
C. Guarnaccia 5 , M. Guida 5 , D. Hatzifotiadou 3 , D.W. Kim 6 , J.S. Kim 7 , S. Kiselev 1 , G. Laurenti 3 , K. Lee 6 , 
S.C. Lee 6 , Ye. Lyublev 1 , M.L. Luvisetto 3 , D. Mal'kevich 1 , A. Margotti 3 , A. Martemiyanov 1 , K. Mikhin 1 , R. Nania 3 , 
F. Noferini 2 - 3 , A. Pesci 3 , M. Ryabinin 1 , E. Scapparone 3 , G. Scioli 2 ' 3 , A. Selivanov 1 , S. Sellitto 5 , R. Silvestri 5 , 

A. Smirnitskiy 1 , Y. Sun 7 , G. Valenti 3 , I. Vetlitskiy 1 , K. Voloshin 1 , L. Vorobiev 1 , M.C.S. Williams 3 , D. Yakorev 6 , 

B. Zagreev 1 , C. Zampolli 2 ' 3 , A. Zichichi 2 - 3 - 4 



> 

o 

(N 

o 

O 

a: 



X 



Institute for Theoretical and Experimental Physics, Moscow, Russia 
Dipartimento di Fisica dell'Universita, Bologna, Italy 
Sezione INFN, Bologna, Italy 

Museo Storico della Fisica e Centro Studi e Ricerche "Enrico Fermi" , Roma, Italy 
Dipartimento di Fisica dell'Universita and INFN, Salerno, Italy 
Department of Physics, Kangnung National University, Kangnung, South Korea 
World Laboratory, Lausanne, Switzerland 

the date of receipt and acceptance should be inserted later 

Abstract. The <f> — > K + K~ decay channel in Pb-Pb collisions at LHC is studied through a full simulation of 
the ALICE detector. The study focuses on possible signatures in this channel of quark-gluon plasma (QGP) 
formation. On a basis of 10 6 collisions at high centrality some proposed QGP signatures are clearly visible 
both in K + K~ invariant mass and transverse mass distributions. The high significance of this observation 
appears to reside heavily on the use of the TOF (Time Of Flight) system of ALICE in addition to its 
central tracking detectors. 
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1 Introduction 



High-energy heavy-ion collisions are considered since a 
long time (see for example yQ) golden events as regards the 
creation and detection of a quark-gluon plasma (QGP). At 
SPS a possible evidence for the production of a new state 
of matter, possibly coinciding with QGP, has been found 
j2||SEll5l- The results from RHIC heavy- ion experiment s 
appear to confirm the QGP hypothesis, and seem to give 
evidence of some additional property of the formed plasma 



n- 

The nature and the characteristics of QGP are in large 
respects even theoretically unknown. This suggests that a 
much larger set of measurements is needed to recognize 
what has been really discovered by the cited experiments 
and to reconstruct the characteristics of QGP. Presumably 
only from the convergence of many different experimental 
signatures can the QGP be fully understood and the QGP 
interpretation of the data established beyond any doubt. 

One of the expected and already detected signatures is 
tied to the s-quark enhancement accompanying (partial) 
chiral symmetry restoration. This should raise the strange 



and multiply strange particle content of the event [7,8,9 



In this respect the <j> particle is of particular interest 
because it is the lightest of the vector mesons with hidden 
flavour. The production and decay rates of these particles 
are strongly influenced by the OZI (Okubo, Zweig, Iizuka) 
rule, that requires the transitions to be described by con- 
nected quark diagrams |15II16| . A first consequence is that 
in case of QGP a strong enhancement of the <j> rate is ex- 
pected as the system has abundance of s-quarks in the 
initial state to feed <j> production. A second consequence 
is the small cross section of the <j) for scattering with non 
strange hadrons. This implies that <p particles are only 
slightly affected by rescatterings with nuclear matter dur- 
ing the expanding phase, rescatterings that could wash 
out the spectral information of the plasma. 

In previous works the capabilities of the future AL- 
ICE experiment at LHC ^TJ in extracting the "normal" 
4> meson signal were discussed, using a full detector Monte 
Carlo simulation 1 8 , IS] ■ The aim of this paper is to present 
the results of a first study, again through full simulation 
of the apparatus, of the potentialities of ALICE in the 
detection of QGP signatures in the </> signal. 
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The study of the <j> in the — > K + K~ channel has 
sonic additional problems with respect to the — > ^ + ^~ 
channel (£ stands for leptons). These are mainly related to 
the fact that kaons are subject to strong interactions, with 
consequent possible rescatterings with nuclear matter, and 
that the <j> —* K + K~ threshold (i.e. twice the kaon mass) 
and the nominal 4> mass are close to each other, with phase 
space problems in the study of possible modifications of 
the mass distribution. 

On the other hand the (j> — > K + K~ decay has a branch- 
ing ratio of 0.49 whereas the value for any leptonic decay 
mode is of the order of 10 -4 . In addition kaon identifica- 
tion is performed in ALICE with low contamination and 
high efficiency over a large solid angle within a wide en- 
ergy range. These circumstances suggest that at least in 
the initial stage of ALICE data taking, the <f> — > K + K~ 
should be the most interesting channel for <j> studies. 

The present work investigates the sensitivity of AL- 
ICE to QGP signatures connected to the <fi decay channel 
into kaons, in the hypothesis that some kind of modifi- 
cations are present in the <j) mass distribution and/or in 
the 4> momentum spectrum in the plasma phase and that 
they survive, even if largely distorted, after <fr and kaon 
propagation. While performing this, particular attention 
is paid to the gain in the <f> — > K + K~ analysis coming 
from the addition to the TPC (Time Projection Chamber 
20 ) of the TOF (Time Of Flight [21]) system for particle 
identification (PID). This in view of the fact that the pres- 
ence of a powerful TOF barrel detector, extending PID at 
medium/high momenta |21U22| . will be a key of ALICE. 

2 Physics hints 

An enhancement of s-quark production in connection with 
QGP formation has been predicted since a long time 
[S]. An enhancement of the <p/uj ratio, as remarked in |23| . 
should point to this effect in a particularly clean way, but 
also an enhancement in the 0/hadron ratio is in general 
expected. 

The s-quark density is readily calculable in case of 
non-interacting QGP, assuming that the equilibrium level 
is reached [2 25 . During the collision however, even if 
QGP is formed, ordinary nuclear matter is also present 
and, in any case, QGP evolves to ordinary nuclear matter 
at freeze-out. In addition an s-quark enhancement can be 
expected also without QGP formation, as a consequence 
of pure nuclear effects, so that only the details of the en- 
hancement as a function of nucleon-nucleon center-of-mass 
energy, ^/snn, could permit to discriminate among QGP 
and non-QGP scenarios. 

It has been suggested that the dependence on ^/snn 
of the s-quark content in case of QGP formation could 
be non monotonic f 26I271IZH ] , this being detectable for ex- 
ample in the s/ir ratio. Experimental confirmation for this 
prediction can perhaps be seen 29,30 from the K/ir mea- 
surements obtained at AGS/SIS [H], SPS [23 and RHIC 

muni. 

Then even though the s-quark enhancement is a very 
interesting item, there are some ambiguities in the evalu- 



ation of the expected amount in case of QGP formation 
and in the interpretation of the enhancement already de- 
tected so far. It is generally believed that this item alone is 
not (and presumably will not be, in the future) sufficient 
to constrain to QGP interpretation; other items must be 
considered in addition. 

In the present work we are concerned with signatures 
of QGP related to the <f> signal alone. This means that we 
will try to estimate the sensitivity with respect to possible 
changes in some properties of <f> mesons, for example in 
their mass distribution, assuming a reasonable scenario 
for s-quark enhancement at LHC energies. 

The 4> signal in the tj> — ► K + K~ decay channel is ex- 
tracted from the large K + K~ combinatorial background. 
Also this background is influenced by the amount of s- 
quark enhancement. However if the K and <p yields grow 
in the same proportion, the <j> signal significance in first 
approximation does not change. A key element in the sim- 
ulations is therefore the 4>/K ratio. 

According to experimental results 12,13,14,29,30 , the- 
re is no sign of a significant increase in the K + /ir + ratio 
going from SPS to RHIC energies, with a kind of plateau 
already starting at the lowest SPS energies. Different is 
the situation for the K~ /n~ ratio, for which a rise from 
SPS to RHIC energies is visible; at the highest energies the 
two ratios, for positively and negatively charged particles, 
are almost equal. 

As regards the <$> meson, a monotonic increase in the 
4>/h~ ratio with ^snn is present [HOj (h~ stands for neg- 
ative hadrons), starting from few GeV energies [HI an d 
going through SPS |3l25l3n| to RHIC energies j3lTl3^] . The 
<j>/K~ ratios found at STAR |32I33I33] and PHENIX [H] 
appear however to be roughly the same as at top SPS 
energies [H]- 

The K/ir ratio within the HIJING generator [12] is 
0.11 aty / T/vw = 5-5 TeV, i.e. a bit lower than the RHIC 
value (~ 0.15) |12II13II14| . The <p/K~ ratio, as measured 
at RHIC 14 , 32 , 33 , 34 , is consistent within errors with the 
HIJING value of 0.16 at 5.5 TeV. 

In our simulations we have used the HIJING values at 
^/snn = 5.5 TeV both for kaon and <j> multiplicities. As we 
are referring to LHC energies, this appears to be a rather 
conservative assumption. In fact the quoted experimental 
results suggest that above RHIC energies kaon production 
presumably saturates, whereas there is no indication in 
this sense for cf> mesons. As a consequence the 4>/K ratio at 
LHC energies is likely to be larger than at RHIC energies. 

During heavy-ion collisions at LHC energies, <f> gener- 
ation and decay happen in a heterogeneous and variable 
medium for which two components can be envisaged. One 
is purely nuclear, consisting of nucleons/hadrons at vari- 
ous number densities and temperatures, but far from T c , 
the critical temperature for the onset of QGP transition; 
the other one is the foreseen QGP component. The latter 
may also include a mixed-phase configuration (occurring 
during plasma expansion) in which some plasma effects 
still persist. To the former component belong in particu- 
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lar all nucleons/hadrons present after chemical freeze-out 1 
is reached. Let us call it in the following the freeze-out 
component. 

To estimate the relative yield of <j> meson decays from 
the two components is a quite complicated task. It de- 
pends on parameters theoretically largely undetermined 
and has no hope of experimental assessment at present. 
In case of a first order phase transition, the estimated du- 
ration of the QGP phase is relatively large (~ 10 fm/c, 
comparable with a vacuum lifetime of <fr mesons of ~ 45 
fm/c) and it can be expected that the two yields are com- 
parable . This should hold true even if the phase transi- 
tion is instead smooth [21], since lattice calculations show 
that in this case the phase transition is a crossover very 
close to a first order one |57) . 

At LHC energies the cf> mesons if decaying inside the 
freeze-out component should have their free-space Breit- 
Wigner mass distribution centred at 1019.5 MeV/c 2 with a 
width of 4.3 MeV/c 2 . Modifications of the resonance mass, 
width and shape can indeed be present in nuclear matter 
at high density, coming from phase space distortions and 
dynamical interactions with the medium |38| . They are 
expected to be observable, but small. 

For decays inside the QGP component the foreseen sit- 
uation is very different. First of all, as a consequence of 
partial chiral symmetry restoration, the mass distribu- 
tion may be at sizeably lower values |39U40II41II42) . The 
exact value of the estimated shift is subject to large un- 
certainties, depending on the assumptions made in the 
calculations. The plasma <j> mass could be 10-20% lower 
than the nominal mass value, but also very near to it. 

Since the kaon mass is expected to change very smoothly 
with temperature 03], the shifted cf> mass can also be well 
below the K + K~ decay threshold in the plasma. This im- 
plies that in principle the width of the resonance should 
decrease in the plasma. However the width is also influ- 
enced by additional effects related to the presence of a 
medium. These could result from an attractive kaon po- 
tential favouring <f) decays |44j (negligible when the <fr mass 
is well below the two-kaon mass) or from <j) interactions 
with surrounding hadrons, such as 7r, K, p and cf> itself, en- 
abling <p decay modes otherwise not present in free-space 
|43II45| and changing the 4> spectral function. As a net 
result an increase of the width up to ~ 10 MeV 43^ can 
fairly be expected. Further few MeV are expected to come 
through the interactions with partons in the plasma |46| . 
The calculated width depends on the effective interaction 
model chosen to describe the plasma. Some choices can 
give even larger widths, of the order of a few tens of MeV, 
as in the case of |47j . 

On top of this, as we are studying <f> decays in the 
<t> — > K + K~ channel, we must consider the effects on the 
observed width of kaon collisions with light mesons, both 
for the freeze-out and QGP components. In- medium kaon 
collisions change the kaon spectral function (off-shell kaon 
final states) and, as a consequence, the <j) decay rate. In 
addition, in-medium (QGP or freeze-out) rescatterings of 

1 When, in the evolution of the hadronic system, inelastic 
collisions cease and particle ratios become fixed. 



on-shell kaons should have an influence on the measured 
width. In principle these rescatterings can even completely 
wash out the signal. In the present work we assume that 
a sizeable fraction of <j) — > K + K~ decays occurring in the 
QGP component is "useful" to produce a signal with an 
additional width, due to off-shell kaon effects coupled with 
kaon rescattering effects, comparable with the original <p 
width. Also in the case of decays in the freeze-out compo- 
nent, a width larger than the normal one has been used 
as a way to describe in-medium effects. 

According to theoretical expectations jSE] as well as to 
experimental hints [T4"ll31ll33| . the transverse momentum 
(or transverse mass) distribution of cj> mesons generated 
inside the freeze-out component at LHC should be harder 
than the distribution found at RHIC, that is with an in- 
verse slope (i.e. temperature) parameter T larger than the 
RHIC measured value: T = 379 ± 50 ± 45 MeV [2], or 
T = 363 ± 8 (stat.) MeV (at 0-5% centrality, with similar 
values for other centralities) or T = 366 ± 11 ± 18 
MeV (minimum bias) |14j . 

On the contrary a much softer distribution is expected 
for </> particles generated inside the QGP component, with 
T ~ T c (independently of the initial temperature 36 ) 
and typical expected values for T c in the range 160-180 
MeV. Due to the small probability of rescattering of the 
<j> mesons after chemical freeze-out [53], the <f> particles 
generated in the plasma should retain their energy distri- 
bution in the plasma even when they subsequently decay 
in the freeze-out component. All this suggests that trans- 
verse mass 2 (mt) distributions are another powerful tool, 
in addition to invariant mass (rrii nv ) distributions, to ex- 
tract evidence for QGP formation from the data and, in 
case, to derive T c from slope measurements [3"o] . 

The expected line-shapes of K + K~ invariant and trans- 
verse mass distributions emerging from the described sce- 
nario are then as follows. 

The rriinv distribution can have two peaks, one local- 
ized at the normal (free space) <j) mass value and the other 
at lower values. The normal peak should correspond to </> 
particles either produced and decaying in the freeze-out 
component or produced in the QGP component but de- 
caying afterwards in the freeze-out one. The lower mass 
peak should be formed by <f> mesons produced and decay- 
ing in the plasma. 

A "bridge" between the two peaks can be expected, 
reflecting the possible mixed-phase evolution from plasma 
to chemical freeze-out. 

The lower mass peak can also be missing; of course 
this may occur when the <fi mass in the plasma is lower 
than the two-kaon threshold. On the other hand, the QGP 
peak can also be so near to the normal freeze-out peak 
to be practically undistinguishable. However, even in this 
case, a signature of QGP formation could be searched for 
in the rrit distribution as we will see later on. A double 
component could/should in fact be present, one (QGP) 
mostly present at low mt values and the other one (freeze- 
out) surely dominant at large values. 

2 ( 1 2\^/^ 

mt = [m inv +pt ) 
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3 Simulation strategy 

Orders of 10 5 -10 6 events turned out to be needed to sim- 
ulate a significant signal for the present analysis. Direct 
full generation, tracking and reconstruction of such a huge 
number of events, corresponding to very high multiplicity 
final states from central Pb-Pb interactions, was estimated 
a rather severe task and simulations have been instead or- 
ganized in two steps, as follows. 

In the first step, some hundred central 3 HIJING [35] 
events, with ^/sWn = 5.5 TeV, were generated, tracked, 
digitized and reconstructed, making use of detailed simu- 
lations of all ALICE detectors (see |18H19| h In these sim- 
ulations a magnetic field of 0.4 T was used. 

From these simulations the relevant parameters for the 
subsequent second step were extracted, as a function of 
transverse momentum (pt), namely the momentum reso- 
lution (in modulus and direction) and, both for the TPC 
and the TOF, the PID efficiencies [181119) . The momen- 
tum is measured by the TPC and the resolutions used are 
detailed in fifflT§\ (for p t in the range 0.2-9.0 GeV/c, the 
relative pt resolution varies between 1% and 2% and the 
angular resolutions go from 5.0 to 0.3 mrad for both polar 
and azimuthal angles); the TPC PID capability relies on 
dE/dx measurements. 

In the second step 10 6 events were generated, each 
event containing an average of 850 charged kaons. This is 
the number expected from HIJING Pb-Pb central events 
at y/SNN = 5.5 TeV in the pseudorapidity (ry) interval [-1, 
!]• 

The kaons are generated with independent random 
choices of rapidity (y) and pt- The y distribution is a 
gaussian with a = 4 [TSITHg . The kaon p t distribution 
is obtained through the mt-scaling method starting from 
the corresponding distribution for pions; for the latter, for 
Pt < 0.5 GeV/c, an m t -scaling spectrum has been assumed 
with T = 160 MeV and, for p t > 0.5 GeV/c, the power 
law parametrization reported in |48| has been chosen. This 
gives a kaon mt distribution with an inverse slope param- 
eter varying with m t ; its fitted value being 309 MeV for 
(m t — m) in the range 0-1 GeV/c 2 , and 390 MeV in the 
range 1-2 GeV/c 2 , where m is the nominal kaon mass. 
Both y and mt distributions are the same as used in |18l 

Eg 

In each event an average number of 85 additional char- 
ged kaons (that is 10% more) were generated according 
to the same distributions to take into account in a simple 
but conservative manner |21j contamination effects mainly 
coming from pions. An average number of 35 (j> — > K + K~ 
decays per event was also added, with both kaons inside 
the indicated 77 range. Given the chosen amount of gener- 
ated kaons, this number corresponds to the already quoted 
4>/K~ ratio of HIJING. 

The (f> particles are similarly generated through inde- 
pendent random choices of y and pt and they are allowed 
to decay isotropically in the center-of-mass frame. Their 

3 The impact parameter range used in HIJING is 0-3 fm. 
The resulting average charged particle density at mid-rapidity 
is (dN ch /dy) ~ 6500. 
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Fig. 1. All-decay <p mass distribution for case A (see text); the 
<j> mass inside the plasma is assumed to be below the two-kaon 
threshold; only a small fraction of <j> particles are allowed to 
decay inside the plasma fireball; a sizeable fraction is assumed 
to decay in a mixed environment corresponding to the evolu- 
tion phase from plasma to dense nuclear matter, which gives 
rise to the so-called "bridge" region. 
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Fig. 2. All-decay (f> mass distribution for case B (see text); the 
cj> mass inside the plasma is assumed to be above the two-kaon 
threshold. 



y distribution is the same as reported above for kaons, as 
in 18T9 . Their p t distribution, differently from |18U19| . 
was not obtained through the m t -scaling method used for 
kaons. It required instead a dedicated approach, according 
to the purpose of the present work where we have tried to 
study the sensitivity to QGP signatures when the input 
4> mass and p t distributions are actually affected by QGP 
formation. To this aim <j> mass distributions different from 
the PYTHIA 0Hj one have been chosen and particular m t 
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distributions, for (j> mesons generated in both the QGP 
and the freeze-out components, have been used. 
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Fig. 3. Invariant mass distribution for all opposite-sign iden- 
tified kaon pairs in 10 6 events. The normalized (see text) back- 
ground distribution relative to equal-sign identified kaon pairs 
is superimposed. The results with/without TOF are compared. 



For each event all kaons were finally processed ac- 
cording to momentum resolutions and PID efficiencies ex- 
tracted from the first step simulation |18U19| . 

Concerning the generation of the </> invariant mass spec- 
trum, two very different situations may appear: one (call 
it case A) in which the cf> mass inside the plasma (referred 
to as QGP c/> mass in the following) is below the two-kaon 
mass threshold and another one (case B) in which it is 
above. 

For case A, the results turn out to be quite insensitive 
to the exact value of the QGP <fi mass as well as to the 
relative importance between the pure QGP part and the 
freeze-out part of the <fi mass distribution. What is cru- 
cial instead is the contribution assumed for the "bridge" 
region, corresponding to the evolution from QGP to chem- 
ical freeze-out. 

In Fig. ^ an example of generated (f> mass distribu- 
tion representative of case A is reported. Here all cf> decay 
modes are included (with the <f> produced in the y interval 
[-1, 1]) and rrii nv is the invariant mass of any decay prod- 
ucts. Two (QGP and freeze-out) Breit-Wigner peaks, each 
with an enlarged width T = 10 MeV/c 2 , are assumed. The 
freeze-out peak is centred at the nominal <f> mass value 
(mo). The additional "bridge" region is taken as a uni- 
form (flat) distribution extending from one peak maxi- 
mum to the other. Normalization is such that the number 



of cj> —> K + K~ decays (strictly speaking, the number of 
useful 4> decays, in the sense explained in SectionEJ is ob- 
tained through multiplication by a factor <P(m)/<l r (mo), 
where <P(m) and ^(mg) are the 2-body K + K~ phase 
space factors, respectively, at <fi masses m and mj. This 
gives zero <fr ^ K + K~ actual decays below the two-kaon 
threshold. The peak around the chosen QGP cf) mass (0.9 
GeV/c 2 ) is, in this example, depressed. The population 
of <j> decays within the "bridge" region is assumed instead 
equal to that of the freeze-out peak, i.e. the peak above 
the KK threshold, at the nominal <f> mass. 

For case B, the results both in the K + K~ invariant 
mass line-shape and in the K + K~ transverse mass spec- 
trum depend on the exact location of the chosen QGP <fi 
mass. 

In Fig. 2 (analogous to Fig. 1) an example of <fi mass 
distribution generated for case B is reported. In this ex- 
ample, the QGP tfi mass peak is at m = 1004 MeV/c 2 , 
which is approximately half way between the KK thresh- 
old and the nominal </> mass. Here the two Brcit-Wigner 
peaks are assumed to be equally populated and the addi- 
tional "bridge" region strongly depressed. 

Turning now to the corresponding m t spectra, these 
are described by fixed inverse slope parameters. For <fi 
particles generated in the freeze-out component, the value 
T = 480 MeV (with a spectrum close to that in [HIIH|) 
has been chosen. This value is higher than the measured 
RHIC values (131^1133) quoted in Sectional thus allowing 
to take into account the effects of a larger initial tem- 
perature at LHC energies. The critical temperature for 
the transition to QGP has been assumed to be T c = 180 
MeV, hence this value was chosen for <f> mesons generated 
in the plasma |36| . 

For each mass distribution used in the present exer- 
cise (Fig. H or EJl, the corresponding m t distribution is the 
sum of only two components. The first, with inverse slope 
T , corresponds to 4> meson decays belonging to the freeze- 
out Breit-Wigner peak; the second, with inverse slope T c , 
to all the other </> decays belonging to the QGP peak and 
the additional "bridge" (if any). As said in Sectional the 
freeze-out Breit-Wigner peak should in principle collect all 
4> decays occurring in the freeze-out, no matter where the 
4> mesons are actually produced (although most of them 
would likely be produced in the freeze-out). Therefore as- 
suming for the whole freeze-out peak a unique tempera- 
ture T corresponds to neglect, as a first approximation, 
the cases in which the (f> is generated in the plasma (or in 
the "bridge") and decays in the freeze-out. 



4 Results and discussion 

Figure [3] shows the inclusive rrii nv distribution for all iden- 
tified K + K~ pairs relative to the 10 6 simulated events 
obtained as described in Section [31 using as input the 
single-peak cf> mass spectrum of Fig. ^ (case A). The AL- 
ICE detector effects are taken into account and two con- 
figurations are compared: one in which the TOF system 
is not present in the detection and PID chain ("without 
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TOF"), the other in which the TOF system is also con- 
sidered ("with TOF"). The difference in the K+K~ pair 
population for the two cases is evident; there is roughly a 
ratio 1.8 in the regions just below and above the nominal 
4> mass. This difference is obviously determined by the en- 
largement, thanks to the TOF, of the momentum region 
where kaons can be identified. 

A significant <f> signal can be extracted after back- 
ground subtraction. For the background we have used 
analogous ui^ nv distributions corresponding to all equal- 
sign identified kaon pairs . The result is shown in Fig. 4. 
This figure actually represents what can happen when the 
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Fig. 4. Extracted signal after background subtraction. The 
assumed <p mass input spectrum is that of Fig. 1 (when the cf> 
mass inside the plasma is below the two-kaon threshold). The 
results with/without TOF are compared. 



QGP 4> mass is below the KK threshold. 

In Fig. 01 again a comparison is made between the 
results with and without the TOF system. The ratio be- 
tween the two peak heights is about 3.3, that is signif- 
icantly larger than the above 1.8 ratio between the two 
background levels (with and without TOF, in Fig.OJ); the 
<f) signal significance is correspondingly enhanced by a fac- 
tor ~ 3.3/V1-8 ~ 2.5. This shows that the larger momen- 
tum acceptance for PID, only achievable when the TOF 
is included in the analysis, reflects itself in a larger sig- 
nal/background ratio, due to the different spectral prop- 
erties of background kaons with respect to kaons coming 
from <p meson decays. 

The fit to the data (with TOF) reported in Fig. H 
is a Breit-Wigner weighted with the appropriate phase 
space factor. The fit parameters, mo — 1018.6 MeV/c 2 
and r = 11.5 MeV/c 2 , are in accordance with the input 
parameters used in this exercise and with the excellent 



mass resolution of the detector |18U19j . A small excess of 
events with respect to the fit is visible, on the left side of 
the peak; this excess comes from </> — > K + K~ decays in 
the "bridge" region defined in Sections [3 and 

Figure [S] shows the K + K~ invariant mass distribution 
for the signal after background subtraction, when the <fi 
mass inside the plasma is above the KK threshold (case 
B of Sect ion 0; this figure corresponds now to the cj> mass 
input spectrum of Fig. 2. Here a double peak structure ap- 




m. 



1.08 

[GeV/c 2 ] 



Fig. 5. Same as Fig. |I] but relative to the 4> mass input spec- 
trum of Fig. 2 (when the <f> mass inside the plasma is above 
the two-kaon threshold). 



pears. One can notice that, with TOF, the QGP <f> mass 
peak is depressed with respect to the freeze-out peak, de- 
spite the fact that the two peaks are equal in the input 
spectrum of Fig. |3 while, without TOF, the two peaks 
remain balanced. This is a consequence of the different 
phase space factors for the two peaks, combined with the 
different PID efficiencies vs. momentum for the TPC and 
TOF detectors. 

When the peaks are too near, they become practically 
^indistinguishable and other ways should be searched for in 
order to disentangle a QGP contribution. Table ^reports 
the <f> signal significance (within ±r/2 from the mean, 
with r = 10 MeV/c 2 ) for both peaks and for various 
choices of the QGP 4> mass. The results are given with 
and without TOF 5 . 



4 The like-sign (K + K + , K K ) background is normalized 
to the K + K~ spectrum in the mass regions below and above 



the (f> peak. 



It can be noticed that the with/without TOF significance 
ratios in this table — from top to bottom — increase from 1.2 
to 1.5 for the QGP peak and decrease from 2.7 to 1.9 for the 
freeze-out peak. Again this is a consequence of the different 
temperatures assumed for the two (QGP and freeze-out) com- 
ponents and of the different performances vs. p t of the two 
(TPC and TOF) detectors. In particular this explains why 
when the QGP and freeze-out peaks are superimposed (last 
table row) the significance gain due to the TOF is lower than 
for the single freeze-out peak of Fig. 0] 
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Table 1. <f> signal significance for case B (see text) in 10 6 
events; the cf> mass inside the plasma is assumed to be above the 
two-kaon threshold and the all-decay <f> mass input spectrum 
to contain equally populated QGP and freeze-out peaks. 



QGP cj> mass 


significance 


significance 


(GeV/c 2 ) 


without TOF 


with TOF 




QGP 


freeze-out 


QGP 


freeze-out 


0.99 


32.8 


41.3 


39.2 


112.2 


1. 


41.9 


37.4 


53.0 


98.7 


1.01 


48.7 


43.3 


74.0 


101.5 


1.02 (nominal) 


67.7 


128.1 



ceptance 6 , just to show how large is the improvement ob- 
tained by adding the TOF. In particular, the region with 
(m t — rriinv) > 0.8 GeV/c 2 cannot be explored without 
the inclusion of this detector |19| . 

In Fig. El one can see that the m t acceptance in the 
configuration without TOF falls off very rapidly for (m t — 
mi nv ) below 1 GeV/c 2 . In these circumstances uncertain- 
ties in mt acceptance would heavily interfere with genuine 
mt spectra and make the QGP discrimination a very diffi- 
cult task, contrary to what happens if the PID capability 
of the TOF is enabled. 

Figure allows to determine the mt region relevant 
for QGP discrimination. The figure refers to the input 



One can see from this table that the signal significance 
on a basis of 10 6 collision events is actually very high. 
In principle, with a statistics of only some 20,000 events, 
a significant 4> signal should already appear. The point 
is however that important handles to disentangle a QGP 
signal are the details of the invariant mass line-shape and, 
as we will see next, of the transverse mass distribution. In 
view of this a superabundant significance of the signal is 
mandatory. 

As said, a promising signature of QGP formation could 
show up in m t distributions. Figure [5] reports the K + K~ 
(mt — rrii nv ) distributions relative to the signal in Fig. ^ 
(referring to the input spectrum of Fig. 1) for the cases 
with and without TOF. The background is subtracted us- 
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Fig. 6. (mt — rriinv) distribution of the signal after background 
subtraction. The assumed (j) mass input spectrum is that of 
Fig. corresponding to a plasma (f> mass below the two-kaon 
threshold. The results with and without TOF are compared. 



ing the same procedure as for Fig. 0] The distributions 
are not corrected for detector K + K~ transverse mass ac- 
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Fig. 7. (nit — rriinv) distribution of the signal after background 
subtraction, in two different regions of invariant mass (see text 
and Fig. QjJ: one around the vacuum value (freeze-out peak) 
and the other one around the plasma value (QGP peak). The 
assumed (j> mass input spectrum is that of Fig.|5| corresponding 
to a plasma <f> mass above the two-kaon threshold. 



spectrum of Fig. [5] and shows the K + K~ (m t — mi nv ) 
distributions (with TOF) for the signal in Fig. [SI after 
background subtraction, in two mi nv regions: one corre- 
sponding to the QGP and the other to the freeze-out Breit- 
Wigner peak. Each region in Fig.UJis selected within ±_T/2 
around the corresponding mi nv peak value. Here again the 
distributions are not corrected for detector K + K~ trans- 
verse mass acceptance. The difference in slope of the two 
distributions is clearly evident up to (m t — mi nv ) ~ 1.5 
GeV/c 2 , with a softer spectrum relative to the QGP peak. 
Above this value the shapes of the distributions become 
quite similar. This indicates that the contamination from 
freeze-out decays within the QGP peak range grows with 
(m t — m^v) and becomes nearly 100% above that value. 



6 This to allow the cleanest possible comparison (i.e. without 
introducing additional biases deriving from the fact that the 
correction for m t acceptance maintains a dependence from the 
input spectra assumed for cj> production). 
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The input spectrum of Fig. 2 has been chosen as an ex- 
ample of a particularly clean situation, to show the princi- 
ple of the discrimination between QGP effects and freeze- 
out effects. Nevertheless, using the present 10 6 event data 
sample and taking advantage of the TOF PID, it can be 
shown that a similar discrimination in m t spectra can be 
achieved for every location of the QGP peak above the 
KK threshold (assuming for dehnitencss equal heights for 
the two peaks in the all-decay rrii nv distribution). When 
the two peaks are very near or exactly superimposed, the 
two mt spectra refer to the same rrii nv regions around the 
nominal <fi mass value; in this case QGP effects are still 
visible as a change in the slope of the spectrum, becoming 
harder above (m t — m inv ) ~ 0.5 GeV/c 2 . 

Such an m t discrimination thanks to the TOF is in 
principle still viable when the QGP <j> mass is below the 
KK threshold, as in the situation of Fig.s 1 and 4, even 
though with smaller statistical significance. The two m„„ 
regions to be compared in this case are the freeze-out peak 
and a narrow region just above the KK threshold. 

As also ordinary nuclear matter can give some shift 
in the resonance mass as well as a distortion in its shape 
|38| . the presence of anomalies in invariant mass distribu- 
tions of the kind described in Figs. 0] and [5] could fail to 
be a compelling indication of QGP formation. Then the 
use of the additional m t information turns out crucial to 
obtain such an indication, whichever the invariant mass 
line-shape is. 

In the results shown above the signal has been ex- 
tracted, as we have seen, through a procedure of uncor- 
rected background subtraction. Background is however 
expected to be uncorrelated only at first approximation. 
As there is no knowledge of particles decaying to K + K~ 
with mass near the nominal (f> mass, the most interesting 
sources of correlated background potentially relevant for 
the present analysis are particles with irp, irK or tttt de- 
cays where one or both decay products are misidentified 
as kaons. 

The probability that the PID procedure misidentifies a 
pion or a proton as a kaon does not exceed the few % level 
|18U19| . This rules out any significant contribution from 
7T7T or irp decays of abundantly produced nearby particles 
such as Kg or A (moreover, for the former, kinematics 
pushes this contribution mainly above the m inv region of 
interest). Hence in practice only ttK decays can have a rel- 
evance. For this kind of decays the best candidate is K°*. 
But again, in addition to a ~ 10~ 2 reduction factor from 
PID efficiency, almost all contamination decays from this 
meson are shifted to mi nv values larger than the nominal 
4> mass due to kinematics. The assumption of uncorrelated 
background appears then fairly safe. 

A last comment regards final-state kaon interactions. 
All the results obtained reside on the assumption that a 
non negligible number of <f> decays have the kaons coming 
out from the plasma with momenta resembling their ini- 
tial momenta (the so called useful <fi decays, see Section 
[2J). To change this number, i.e. to change the relative im- 
portance of the QGP and freeze-out components in the 
all-decay input spectra (Figs. fl] and |2J , implies to change 



the significance of QGP discrimination. With a reduced 
number of useful decays this discrimination can of course 
become impossibile. 

On the other side it seems quite difficult to reproduce 
the described signatures of QGP formation through re- 
interactions of kaons (or of pions and protons misidenti- 
fied as kaons) coming from decays of resonances other than 
the <p. Re-interactions will presumably spread kaon (and 
other particle) momenta and distort their spectra, but it 
seems unlikely that they would be able to mimic combined 
signatures of the kind described herein in both the invari- 
ant mass and transverse mass distributions. From this one 
can infer that if such signatures are observed they would 
strongly point to an origin not easily reducible to final- 
state interaction effects, and contribute to compel to a 
QGP scenario. 



5 Conclusions 

In this work the sensitivity of ALICE to QGP signatures in 
Pb-Pb collisions has been investigated through full detec- 
tor simulation by studying the <fi — > K + K~ decay channel. 

The impact of QGP formation on two quantities, the 
mass and the transverse momentum of </> mesons produced 
inside the plasma, has been considered. 

On a basis of 10 6 central events, a discrimination of 
QGP effects has been shown to be possible for various 
choices of <f> mass inside the plasma, even below the KK 
threshold. This can be achieved by looking at secondary 
peaks or distortions of the invariant mass distribution of 
the extracted signal in the region below the nominal <f> 
mass and at the spectral properties (transverse mass dis- 
tribution) of this signal. 

The combined study of both the rrii nv and the m t dis- 
tributions of <fi mesons has been identified as a powerful 
tool of investigation. 

The significance of the results depends on various the- 
oretical unknowns such as the QGP <fi mass value, the du- 
ration of the QGP phase, the evolution mechanism of the 
system from plasma to chemical freeze-out, the particle 
ratios and spectra, among others. 

Experimentally this significance would strongly benefit 
from the TOF system for particle identification. The over- 
all </> signal significance is actually found to be enhanced 
by a factor 2-3 thanks to the TOF; moreover, this gain 
in significance strongly increases with m t while, in the 
configuration without the TOF, no sensitivity is present 
when (m t — rrii nv ) is above 0.8 GeV/c 2 . This means that 
the effectiveness in QGP discrimination in terms of mt 
appears to depend crucially on the presence of the TOF 
PID capability. 
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